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We study Klein tunneling in polariton graphene. We show that the photonic spin-orbit coupling 
associated with the energy splitting between TE and TM photonic modes can be described as an 
emergent gauge field. It suppresses the Klein tunnelling in small energy range close to the Dirac 
points. Thanks to polariton spin-anisotropic interactions, polarized optical pumping allows to create 
potential barriers acting on a single polariton spin. We show that the resulting spin-dependent Klein 
tunneling can be used to create a perfectly transmitting polarization rotator operating at microscopic 
scale. 
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Emergent physicJ^ in solid-state systems with effec¬ 
tive Hamiltonians mimicking the behavior of less acces¬ 
sible systems has become a very productive field of re¬ 
search. The work of SemenoflP was a starting point for 
the research on the analogs of effective electrodynamics 
in solids, which has culminated with beautiful works on 
graphen^, whose practical properties, such as the car¬ 
rier mobility, are very much affected by the effect pre¬ 
dicted in high-energy physics, but inaccessible there - 
Klein tunneling, leading to perfect transmission through 
a potential barrier by particle-antiparticle conversion. 

Although graphene is a very promising, rich in 
effect^^, and popular system, offer ing wide possibili¬ 
ties for the study of emergent physics^, it has its own 
limitations: many measurements are indirect, its struc¬ 
ture is fixed and its parameters (such as the band gap) 
are difficult to tun^. This is where different types 
of artificial graphen^^ come into play. They ca n be 
based on di fferen t particles: atoms, both fermionicfHE^ 
and bosonicP^^, or on photons and photonic quasi- 
particle^I^®^, confined in a 2D honeycomb potential. 
At the lowest level of approximation, the single parti¬ 
cle Hamiltonian is similar with the one of graphene, and 
it typically results in the same type of dispersion, charac¬ 
terized by the presence of the famous Dirac cones. These 
different systems offer a very wide tunability. The lat¬ 
tice parameters can be modified, spe cific types of spin- 
orbit interaction can be implementecP^^^. The physics 
of both bosonic and fermionic interacting systems can 
be addressed. In photonic systems, time-dependent per¬ 
turbations can be used to create Floquet topological 
insulatorJ^. Photonic systems in general allow a unique 
direct access to the time and spatial evolution of the wave 
functions (for example, the Bloch functions of a lattice in 
real a nd rec iprocal space), with very simple experimental 
mean^^^^. 

Recently, a 2D honeycomb lattice has been im¬ 
plemented for interacting photons, the exciton- 
polaritons (polaritons]P^. These quasiparticles appear in 
microcavitie^^ in the strong coupling regime between 
the quantum well excitons and the cavity photonJ^. 
They combine light effective mass with strong interac¬ 


tions. Their bosonic character provides the possibility 
for Bose-Einstein condensatiorP^, while the two spin 
projections make it possible to describe them within the 
pseudospin formalisnP^ where the pseudospin dynamics 
is described by its coupling to effective magnetic fieldJ^. 
The spin-anisotropic character of the inter act ionPSl 
together with various controllable effective fields offer 
a large variety of spin (polarization) effects for spin- 
optronics. The fabricated polariton graphene is based 
on a lattice of coupled micropillar^. This system is 
characterized by a spin-orbit interaction (SOI) acting 
on the real polarization of the photonic eigenstates. 
This SOI is induced by the energy splitting between the 
TE and TM polarized eigenmodes. It makes polariton 
graphene suitable for the realization of the optical 
spin Hall effecfP^. As noticed in different contexts, the 
specific angular dependence of the TE-TM induced 
SOI induces chirality, whic h can generate stationary 
photonic spin currentJ^^^^. When combined with a 
Zeeman effective field in polariton graphene, it leads to 
the for mation of a polaritonic analog of a Z topological 
insulato j^^ l ^^ l 

In this work, we propose and analyze an experimental 
scheme to study Klein tunneling in polariton graphene. 
The scheme of the proposed experiment is shown on 
Fig. 1(a). The potential barrier is created by a non¬ 
resonant optical pumping (blue), which populates an ex- 
citonic reservoir. This excitonic reservoir creates via the 
exciton-exciton interactions a potential barrieJ^ affect¬ 
ing polaritons. The useful signal consists of polaritons 
optically injected close to the Dirac point of the disper¬ 
sion by quasi-resonant pumping (red). The height of the 
barrier can be very precisely, controlled as demonstrated 
experimentally^^H^. Our study will include the role of 
the specific SOI present in the system, and also the free¬ 
dom offered by the spin-anisotropic interactions between 
polaritons, which allows to optically create a potential 
barrier acting only on one specific spin componenfPSi. 
Klein tunneling in the presence of Rashba like spin-orbit 
coupling has already been studied both for electronJ^ 
and Bose-Einstein condensateP^. Qualitatively similarly 
with their findings, we show that the TE-TM induced 
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Figure 1: (color online) a) Scheme of the experiment. Pat¬ 
terned microcavity shown in grey, pump 1 (red) creates prop¬ 
agating polariton, pump 2 (blue) controls the barrier; b) 
Particle-antiparticle conversion in Klein tunneling: the dis¬ 
persion is shown in red, the potential in blue, and the parti¬ 
cle in green; c) Scheme of the Klein tunneling configuration. 
The indices z,r,t stand for incident, reflected and transmitted 
waves. The barrier is shown by blue hatching. 


SOI considerably suppress the Klein tunneling in a nar¬ 
row energy range around the Dirac point. We show 
that this field can be represented as an emergent gauge 
field. In the second part of the manuscript we propose 
a scheme of ” Klein polarization rotator” which allows to 
optically control the polarization emitted by a the po¬ 
lariton graphene: if a potential barrier is created only 
for one of the two spin components, it affects the phase 
of this component, and thus the orientation of the linear 
polarization, while the transmission remains perfect due 
to the Klein tunneling. 


I. TIGHT-BINDING DESCRIPTION OF 
POLARITON GRAPHENE 


Before addressing Klein tunneling, we need to discuss 
the system, which provides the effective Hamiltonian re¬ 
quired to achieve this phenomenon. Polariton graphene 
inherits the most important properties of real graphene. 
The simplest tight-binding model of electronic states in 
graphene^ neglects the spin degree of freedom. Tak¬ 
ing into account only the nearest-neighbor tunneling de¬ 
scribed by a constant J, the Hamiltonian in the basis of 
the two atoms A and B of the unit cell forming the lattice 
writes as follows: 


Hk = 




( 1 ) 


where /k = exp(—ikd^^.) is a sum over the 3 near¬ 

est atoms. This model gives rise to an emergent Dirac 
equation for electrons at the special points of the dis¬ 
persion (called K and K’ or simply Dirac points) lo¬ 
cated at the corners of the Brillouin zone. In this equa¬ 
tion, the emergent ’’spin” 1/2 corresponds, in fact, to 
the sublattice degree of freedom - atoms A and B (the 
coordinates in the reciprocal space are usually modified 
for convenience with respect to the full Hamiltonian as 
kx ^ ky ^ky y kx^ - 


= = ( 2 ) 

where vp = ^Ja/2h is the Fermi velocity, replacing the 
light speed in the original Dirac equation. For a bosonic 
particle, it should be considered simply as a parameter 
of the dispersion. The index of the dispersion branch is 
defined as a = sign{Ekin) = il, Ekin being the kinetic 
energy with its zero given by the energy position of the 
Dirac points. By the analogy with the Dirac equation, 
excitations with the branch index ±1 are called ’’parti¬ 
cles” and ” anti-particles” respectively. In order to apply 
this Hamiltonian for exciton-polaritons, formed by strong 
coupling of two dimensional excitons and photons placed 
in a honeycomb potential, we need to make several ap¬ 
proximations. The transverse (lateral) dynamics of pho¬ 
tons in a cavity can be described using the Schrodinger 
equation. Restricting the consideration only to the lower 
polariton branch and using the parabolic approximation, 
valid close to the bottom of this branch, it becomes pos¬ 
sible to apply the tight-binding approach, provided that 
the lattice sites are not too small (to avoid large wavevec- 
tors) and the confinement is strong enough, so that the 
bandwidth is smaller than the energy difference between 
the two first energy states of an isolated lattice site. 

Periodic potential acting on photons and polaritons 
out of planar microcavities can be realized in different 
waysSlHlI^ the main example being however based on a 
lattice of micro-pillars obtained by patterning of a planar 
microcavity for which the above mentioned approxima¬ 
tion is well realized. 


A. Spin-orbit coupling in polariton graphene 

Since we are going to deal with the polarization degree 
of freedom, we need to rewrite the tight-binding formal¬ 
ism, accounting for the two spin projections of polaritons, 
corresponding to right- and left-circular polarized pho¬ 
tons. For this, we will have to work with bispinors instead 
of spinors. The corresponding tig ht-bin ding Hamiltonian 
for polaritons has been derived iipUMl^ together with ex¬ 
tra terms responsible for the spin-orbit coupling, which 
will be discussed below. In the absence of spin-orbit cou¬ 
pling, the two circular components are completely inde¬ 
pendent, and the Hamiltonian is simply a combination of 
two tight-binding graphene Hamiltonians: 
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- the wavefunctions of the two sublattices and two spin 
components. 

Although we are using the wavefunction description of 
polaritons, based on the Schrodinger equation, which is 
possible due to the quantization of photons in the mi¬ 
crocavity in the growth direction, one should not forget 
that polaritons are formed from photons, described by 
the Maxwell’s equations. It is natural to treat the elec¬ 
tromagnetic waves on the basis of TE and TM eigen- 
modes, which are split in energy in presence of different 
media. For polaritons, the TE-TM splitting has been 
thoroughly discussed since Ref.SSl, and many interesting 
effects based on this split ting h ave been demonstrated 
(e.g. optical spin Hall effecfP^^, or acceleration of emer¬ 
gent magnetic monopole^^. In confined structures, the 
TE-TM splitting can be enhancecP^ with respect to the 
planar cavities. For polariton graphene, it has already 
been the s ubject o f extended experimental and theoret¬ 
ical studieJ^®^^^. Therefore, we need to take it into 
account in our description of Klein tunneling in polari¬ 
ton graphene. The Hamiltonian including the TE-TM 
induced spin-orbit coupling can be written as: 


Hk = 


0 Fk 


A 


0 


( 4 ) 


where the block matrices are defined as: 


Fk 


{ hJ fisj\ 

hj)' 


( 5 ) 


Here, J is the polarization-independent tunneling coef¬ 
ficient, whereas is the SOI-induced polarization de¬ 
pendent term, which can be up to 10% of J. Physically, 
it means that the polariton pseudospin rotates around 
the effective field during the tunneling process. The 
wavevector-dependent complex coefficients /k,/^ are de¬ 
fined by the sum over the nearest neighbors: 

3 3 

/k = E exp(-ikdy^.), /^ = y] exp(-i [kd^^. =F ), 

i=i i=i 


Since the Klein tunneling occurs in the region of the 
reciprocal space close to the Dirac point, we can approx¬ 
imate the expressions for /k,/^, keeping only linear in 
k terms around this point, which leads to the following 

expression: H = hvp {^x^x + ^ {^y^y ~ ^xSx)i 

where A = ZbJj2 is the measure of the spin-orbit cou¬ 
pling, Sx and Sy being the pseudospin operators acting 
on the real polarization of particles. Without the loss of 
generality, we consider only one of the two Dirac points. 
The matrix form of the Hamiltonian in this approxima¬ 
tion can be written as: 


/ 0 0 

0 0 

hvpikx^^ky) 0 

\ 2A hv-p{kx^^ky^ 


hvrikx - iky) 

0 

0 

0 


2A 

hvrikx - iky) 

0 
0 


(6) 


The 4 branches of the dispersion described by this 
Hamiltonian are parabolic for low wavevectors because 
of the spin-orbit coupling, and two of them are split-off 
by A. 

E = ±A± ^/A‘^ + {hvFkf (7) 

This result is the consequence of our approximation, 
valid only at intermediate values of wavevector k relative 
to the Dirac points. In fact, the spin-orbit coupling in its 
full form leads to the trigonal warping of the dispersion. 
However, the typical scale of the trigonal warping is so 
small that it can be safely neglected^. 


B. Emergent non-Abelian gauge field 

Within the approximation of intermediate wavevec¬ 
tors used above for the description of the spin-orbit 
coupling in polariton graphene, the Hamiltonian H = 

Jivp (^^xkx ^yky^ + ^ i^y^y ~ ^x^x) can be reformu¬ 
lated to become mathematically similar to that of a 
charged Dirac particle in presence of a vector poten¬ 
tial of a gauge field. This is possible thanks to the re¬ 
duced symmetry of the effective Dresselhaus spin-orbit 
coupling close to the Dirac point, as compared with the 
TE-TM field in the T poin^^. The Dirac Hamiltonian 
for a charged particle read^^ 











4 


where A is the vector potential of the electromagnetic 
field (B = rot A). As for the Schrodinger equation, the 
momentum p is simply replaced by p — e/cA. Here, e is 
the electron charge and c is the speed of light, which in 
the case of the effective graphene Hamiltonian is replaced 
by the Fermi velocity vf, determined by the coupling J 
in the tight-binding model we consider. 

To have the analogy with the Dirac equation, we need 
to define the effective vector potential as: 


A = - 


As' 

he 


( 9 ) 


where we have inverted one of the axes of the polariza¬ 
tion pseudospin: s'x = — s'y = Sy. In the new coordi¬ 
nates of the polarization pseudospin space, the Dressel- 
haus field is converted into the Rashba field. Using this 
vector potential, the polariton graphene Hamiltonian can 
be written in the gauge field representation: 


where the components of the vector potential and 
Ay do not commu te wit h each other. Thus, the gauge 
field is non-AbeliaiP^^^. Written in this form, the equa¬ 
tion for bispinor neutral particles is reduced to a usual 
Dirac equation for spinor charged particles in an emer¬ 
gent gauge field. Since the vector potential contains the 
polarization pseudospin operators Sx^y^ its spatial distri¬ 
bution depends on the current distribution of the polar¬ 
ization pseudospin, determined by the polarization spinor 
part of the wavefunction. Thus, the emergent electro¬ 
magnetic field texture is defined by the polarization of 
polaritons. This leads to interesting effects, such as the 
lensing by an impenetrable defect, as was shown for other 
configurations with polaritonic emergent gauge fieldJ^. 
Detailed study of the consequences of the emergence of 
the gauge field for Dirac equation is a subject for future 
works. 


II. KLEIN TUNNELING 

The description in terms of the Dirac Hamiltonian 
for excitations in honeycomb latt ices is well established, 
and Klein tunneling in graphen^^^^^, as well as with 
atomic condensates in optical latticeJ^, has already 
been demonstrated experimentally. Although th is p he- 
nomenon is described in many review papers^^M]^ 
mathematical description in the simple case is used as 
a basis for the spinor case considered in this manuscript. 
This is why we briefly revisit the simplest scalar case in 
the subsection below. 


A. Scalar case 

Qualitatively, Klein tunneling for massless particles 
consists in a perfect transmission through potential bar¬ 
riers because of the particle-antiparticle conversion possi¬ 
ble for the Dirac equation: a particle with energy E <V 
turns into an antiparticle with energy E' = E — V^ prop¬ 
agating in the same direction, as shown in Fig. 1(a) by a 
green arrow (dispersion is shown in red). The backscat- 
tering is suppressed because of the pseudospin conserva¬ 
tion: the potential cannot change the lattice pseudospin, 
because it acts identically on both components A and B. 
Of course, this is true only for a particular propagation 
direction or in Born approximation. 

For the general description of Klein tunneling, we need 
first to write the spinor solution of the Dirac equation 
in the A, B atom basis for an arbitrary propagation 
direction 


T = 

V2 


ik.r 


ae 


i(j) 


( 11 ) 


The propagation direction with respect to the horizon¬ 
tal axis is given by the angle (j) with tan(0) = kyjkx as 
shown on the Fig. 1(b). Let us now consider the inci¬ 
dence of a wave defined by Eq.(ll) on a potential barrier 
of a height V located at x = 0 (uniform along y). 

The energy of the particle given hy E = hvF\k\ in the 
Klein tunneling regime is smaller than the barrier height 
V. The incident, reflected, and transmitted wavefunc- 
tions can be written as: 



The continuity of the wavefunction in x = 0 imposes the 
following constrains on the two components of the spinor: 


ae'^y'^e'^ 


+ re'^^yy 
-h rae'^yye'^^ 


te^Kv 

a' te'^yy e'^^ 


(13) 


The invariance in the y direction imposes the conserva¬ 
tion of the wave vector ky = ky = ky. Since E < V, a = 1 
and a' = —1. This allows to determine the remaining un¬ 
knowns: (l)r = TT — (pi and —{E — Vq) sin(0t) = Esin(0^). 
The equations become: 


{ 



Which leads to: 


t 


(14) 


g ^ 


2 cos{pi) 

g — g*0t 


(15) 
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The expressions become particularly simple for V = 
2E, giving Esin((/)t) = Esin((/)i), or simply = 

IT — (pi. The expressions for the reflection and tunneling 
amplitudes and intensities become: 


r ^ _ g2((/)i+7r/2) sin(0^) , t = cos{(pi) 

\ i? = sin^(^i) , T = cos'^{4>i) 

(16) 

One can see that indeed, for normal incidence (pi = 0 
the reflection is suppressed, and the transmission is T = 
1. This is the famous Klein tunneling effect, relying on 
the particular shape of the dispersion. 


B. Klein tunneling in the spinor case 

The extra degree of freedom given by the polarization 
of light allows to define a second pseudospin, indepen¬ 
dent from the pseudospin associated with the lattice sites 
A and B. This second pseudospin corresponds to the 
Stockes vector of light, and it is related to the compo¬ 
nents of the spinor in the circular basis as: 

5 , = 

Sy = (17) 

5 , = (nf-rf^)/2 

In this section, we are going to consider the specific 
case of a potential barrier present only for one of the two 
spin components (say, cr+), and absent for the other spin 
component (cr“), which can be realized because of the 
spin-anisotropic interactions of polaritons and by using 
circularly polarized resonant or non-resonant pumping 
to create the potential barrieJ^ (see Annex I). Since the 
two wavefunctions are uncoupled (we neglect the spin- 
orbit coupling in this section), the Klein tunneling (with 
the corresponding modification of the wave function) will 
occur only for one component, whereas the other will be 
just freely propagating. It is easy to write the solution 
of the Dirac equation in this particular case, combining 
the homogeneous solution for (j~ with the solution ex¬ 
hibiting the particle-hole transition for cr+. In the bar¬ 
rier region, the action of the spin-polarized potential is 
qualitatively similar to that of a magnetic field causing a 
Zeeman splitting between the circular components: the 
linear polarization will precess around this field. 


The wavefunction of a linearly polarized state 
propagating in a particular direction reads: ^ = 

(l, 6 *^"^+^^)^, where tan 0 = ky/kx as be¬ 

fore, and 0 is the relative phase which determines the 
orientation of linear polarization (or the direction of the 
pseudospin), while a gives the sign of energy (particles or 
holes). Since the two polarization components are essen¬ 
tially independent in the absence of spin-orbit coupling, 
it is useful to write the corresponding spinors separately: 
T+ = (T;^,T^)^, = (T;^,T^)^. The incident, re¬ 

flected and transmitted wave functions for both polariza¬ 
tion components are given by the following expressions, 
based on the previous results for the reflection and trans¬ 
mission coefficients in the Klein tunneling regime for the 
< 7 + component (still assuming E = V/2 for simplicity; 
fc=|k|): 




1 ^ik{x cos (p+y sin cp) 



sin <p ^i(k(—x cos <p+y sin (p)-\-(p-\-7r/2) 

V2 ® 

0 



cos <p i{k(—x cos (p-\-y sin (p)-\-(p) 

V2 ® 



_ 1 COS (/)+?/sin (/)) 



(18) 


There is no reflected wave for , because the barrier is 
present only for cr+. Moreover, in the case of normal inci¬ 
dence, there is no reflection for cr+ as well, as can be seen 
from the sin(0) factor in T+. In this case, the transmis¬ 
sion for both polarizations is equally perfect, but a rela¬ 
tive phase appears between them, as can be seen compar¬ 
ing and for 0 = 0 , = exp(—2i/cx). Since 

the relative phase between the spin components is what 
determines the orientation of the linear polarization, let 
us calculate it explicitly by analyzing the polarization 
pseudospin in the general case of 0 7 ^ 0. We have: 




(19) 


We obtain 8 pseudospin components : x and y on the 
A and B sub-lattices to the left and to the right of the 
barrier. 
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cL,A 


qL,A 


qL,B 


qL,B 

Oy 


^ (1 — sin 0 sin (0 — 2kx cos (/>)) 
— ^ sin (p cos (0 — 2kx cos p) 

^ (1 — sin (p sin (0 + 2kx cos p)) 
^ sin p cos{p + 2kx cos p)) 



qR,A 


qR,B 

qR,B 

Oy 


I COS p cos{p — 2kx cos p) 

— ^ cos p sm{p — 2kx cos p) 
I cos p cos{p + 2kx cos p) 

I cos p sin (0 + 2kx cos p) 


( 20 ) 



X (arb. units) X (arb. units) 


Figure 2: (color online) Pseudospin components for normal 
0 = 0 (panel a) and oblique 0 = tt/G (panel b) incidence. Lin¬ 
ear polarization always rotates in the right half-plane (barrier 
region). 


One can easily see that these expressions verify the 
continuity at x = 0. From them, one can plot the lin¬ 
ear polarization degree in any basis (horizontal/vertical 
or diagonal/anti-diagonal) for any angle of incidence, 
or simply plot the orientation of the linear polarization 
plane r] = 1/2 arctanS'^/S'^,. We note that the linear po¬ 
larization for the two sublattices A and B is not the same: 
the diagonal polarization can be different for 0 7 ^ 0. The 
main consequence is the existence of spatial variations of 
the linear polarization on both sides of the spin-polarized 
barrier. 

In the case of normal incidence, the expressions are 
strongly simplified: 

r = \ I =S^’^ =icos(2M 

\ 57^ = 57-® =0 \ S'f ^ = 5 ®-® = i sin( 2 fca;) 

( 21 ) 


In this case, spatial oscillations of the linear polarization 
are observed only to the right of the barrier, in the trans¬ 
mitted wave. 


Fig. 2 presents the results of the calculations of the 
pseudospin components according to the equations ( 20 ) 
for normal incidence in panel (a) and for oblique inci¬ 
dence in panel (b). One can see that the polarization 
always rotates in the region of the barrier, whereas in 
the left half-plane the rotation depends on the angle of 
incidence. Panel (a) also demonstrates the full transmis¬ 
sion, signature of the Klein tunneling regime. 


The above description captures the important speci¬ 
ficity of polariton graphene, that is, the presence of 
two different pseudospins, and demonstrates its conse¬ 
quences, but for the moment we were neglecting the spin- 
orbit interaction, which can play an important role for 
Klein tunneling, as shown previously for other systemJ^H. 


C. Suppression of the Klein tunneling in presence 
of spin-orbit coupling 


In this section, we consider the same as previously, 
namely a potential barrier acting only on one spin com¬ 
ponent, but including the photonic spin-orbit coupling. 
The corresponding dispersion obtained in the section 
I, parabolic at small wavevectors and linear at larger 
wavevectors, is shown in Fig. 3(a). The corresponding 
eigenvectors, numbered in the order of increasing energy 
for a fixed wavevector, are as follows: 



hvjpk 


hvFk —id) 


hvpk —id) 


hvpk ^—id> 


= 

El ^ 
hvpk ^i(j) 

, «'2 = 

E2 ^ 
hvpk ^id) 

, 4'3 = 

E3 ^ 
hvpk ^id) 

, «'4 = 

E4 ^ 
hvpk ^id) 

(22) 


El ® 


E, 


Es ® 


E4 



V 1 / 

V 1 / 

V 1 / 

V 1 / 


To find analytically the reflection and transmission co¬ 
efficients we use the same approximations as before: the 


potential barrier is invariant in the Y direction, the en¬ 
ergy is one half of the barrier height E = V/2. Another 








































7 



Figure 3: (Color online) a) Dispersion in the vicinity of the 
Dirac point in presence of spin-orbit coupling, b) Transmis¬ 
sion and reflection coefficients for normal incidence 


is not suppressed, and therefore the results obtained in 
the previous subsection for decoupled spin components 
remain valid. In realistic structures, the energy difference 
between the split-off bands is comparable with the TE- 
TM splitting magnitude at the wave vector of the Dirac 
point. It depends completely of the structure geometry 
but typically ranges between a few tens and 100 //eV. It 
is therefore of the order of the mode linewidth in good 
quality samples. We therefore expect the Klein tunnel¬ 
ing suppression to be an observable but relatively weak 
effect. 


assumption is linked with the choice of the initial energy 
branch, determining the eigenvector. One should take 
into account that transitions between the branches oc¬ 
curring in presence of a potential V should conserve the 
polarization in the vicinity of the barrier edge, which al¬ 
lows only 3^1 and 4^2 transitions. We assume that 
the incident wave corresponds to the 4th branch (posi¬ 
tive and split-off). This gives us the following system of 
equations: 

1 + r = t 

(23) 

hvFk r^id) 

E ^ / E' 

With this system we determine three parameters : r, 
t and (pf 


III. KLEIN POLARIZATION ROTATOR 

The rotation of polarization on a micrometer scale in 
the absence of any backscattering (for normal incidence) 
allows to use the proposed structure as a polarization ro¬ 
tator, which can be called ’’Klein polarization rotator” or 
’’Klein waveplate”. The deterministic control of polari- 
ton spin, associated with the other opportunities offered 
by polariton graphene structures, such as one-way surface 
states, makes of it a promising platform for spin-optronic 
applications. 

To check our analytical predictions, we have performed 
a numerical simulation based on a spinor Schrodinger 
equation for polaritons where the honeycomb confining 
potential U{x^y) is taken into account. We consider 
both the situation with and without spin-orbit coupling. 
Without spin-orbit coupling the equation reads: 


t = 


ke~^^ — ktc^^t 
2kcos(j) 


ke~^^ — kte^^t 

k 

4>t = TT + arcsin( —sin0) 
kt 


(24) 

(25) 

(26) 


In the case of normal incidence, the expressions can be 
greatly simplified. 


k — kt 2k 

k + kt' k + kt 


(27) 


We can thus determine the reflection and transmission 
probabilities R = |rp and T = 1 — R, and analyze how 
they are affected by the spin-orbit coupling A. Since the 
Klein tunneling is associated with the linear dispersion, 
it is natural to expect the suppression of the transmission 
for the parabolic part of the dispersion. The results of the 
calculations are represented in Fig. 3(b) as a function of 
wavevector (plotted ioi h = vp = A = 1,0 = 0). Within 

our assumptions, kt = \/ 4A^ + 4Av^A^ k^ 

At lower wavevectors, the transmission becomes 
strongly suppressed, because the Klein mechanism does 
not protect it anymore. However, for higher wavevec¬ 
tors, where the dispersion is linear, the Klein tunneling 




(28) 


+Poe' 


(■—-n)^ 


^iCkr—ujt) 


where ^(r) = ^^(r),^-(r) are the two circular com¬ 
ponents of the wave function, m is the polariton mass, 
r = 25 ps the lifetime. Since the calculation is per¬ 
formed without the tight-binding approximation, only 
the polarization pseudospin remains. We have taken 
777, = 5 X 10“^mo, where mo is the free electron mass. 
Pq is the amplitude of the pumping (identical for both 
components, corresponding to horizontal polarization), 
the size of the spot a = 5 ym in the X direction and 
40 /im in the Y direction. The result of the simulation 
is shown in Fig. 4, demonstrating the inversion of the 
linear polarization degree just after the barrier (located 
at X = 0), which appears in blue on the figure. 

We see that the numerical simulations confirm the ana¬ 
lytical predictions. Indeed, a barrier for one polarization 
component can be considered as an effective magnetic 
field in the Z direction, because it creates a ’’Zeeman” en¬ 
ergy splitting between circular polarization components. 
This effective field rotates the linear polarization in the 
plane, which is at the origin of the observed effect, while 
the Klein tunneling regime provides the suppression of 
backscattering. The average value of backscattered in¬ 
tensity, mostly due to the finite size of the beam leading 
















Figure 4: (color online) Linear polarization degree in a 2D 
polariton graphene. The pump is located at the left edge of 
the hgure, polaritons propagate to the right. A barrier for cr^ 
component is located at A = 0, which leads to the inversion 
of the polarization degree at A = 10 /xm. Color shows linear 
polarization degree pi = {Ih — Iv)/{Ih + /v)- 


to the violation of the normal incidence condition on the 
edges, did not exceed 1%. 

In order to check that the operation of the ” Klein wave- 
plate” remains possible in spite of the spin-orbit coupling 
existing in all real structures, we have repeated the nu¬ 
merical simulations, but with the new Hamiltonian and 
energy dispersion, adding the spin-orbit couping (with¬ 
out applying the linear approximation required for ana¬ 
lytical calculations). For this, an extra term was added 
into the Hamiltonian of the spinor Schrodinger equation: 

+ (29) 

The TE-TM splittin^^ is described by the parameter 
[3 = li? /4m where rapt are the effective 

masses of TM and TE polarized particles respectively and 
m = 2 {mt — mi) /mtmi. The results of the simulations 
are presented in Eig. 5 (a) and (b), the difference between 
them being the value of Eor Eig. 5(a) we have taken 
the typical value of 5% difference between the masses, 
whereas for Eig. 5(b) this value was multiplied by a factor 
5, corresponding to an artificially high value of TE-TM 
splitting, never achieved in real cavities. 

Once again the numerical simulations confirm the an¬ 
alytical calculations. The typical spin-orbit coupling ex¬ 
isting in polariton graphene is too small to create any 
pronounced effects in the configuration of our simulated 
experiment: the ’’Klein waveplate” on Eig. 5(a) oper¬ 
ates as expected, in spite of the non-zero value of /3 and 
small parabolicity of the branches: a region of inverted 
polarization (blue) appears after the barrier. We did not 
observe any significant increase of the backscattered in¬ 
tensity (1%) with respect to the case without spin-orbit 
coupling (Eig. 4). On the other hand, an artificially en¬ 
hanced value of the spin-orbit coupling (3 leads to the 


20 


E 

3 0 



X(|?m) 


20 -20 


X(|?m) 


Figure 5: (color online) Spatial images calculated by solving 
the full spinor Schrodinger equation for the ” Klein waveplate” 
with TE-TM splitting, showing inear polarization degree with 
a) normal TE-TM splitting, b) TE-TM splitting x5. Color 
shows linear polarization degree (same as fig.3). 


suppression of Klein tunneling (20% of backscattering of 
the circular component with barrier) and the associated 
polarization rotation, as we can see in Eig. 5(b): there 
is no region with pronounced inversion of polarization in 
the right part of the figure. Eurthermore, the linear po¬ 
larization degree is not positive in the left part any more, 
which is a signature of a strong reflection on the barrier, 
together with polarization inversion. 


IV. CONCLUSIONS 

To conclude, we have studied the Klein tunneling in 
polariton graphene, taking into account the spinor prop¬ 
erties of polaritons, including the spin-anisotropic inter¬ 
actions and the spin-orbit coupling. We have shown that 
while the interactions allow to exploit the Klein tunnel¬ 
ing effect for the creation of a micron size polarization 
rotator without reflection on its surfaces, the spin-orbit 
coupling might perturb its operation. However, the effect 
of this coupling for realistic structures remains negligible. 

We acknowledge the support of ANR GANEX 
(207681) and Quandyde (ANR-ll-BSlO-001), and EU 
ITN INDEX (289968) projects. 


V. ANNEX I 

Polaritons, due to their mixed light-matter nature, ex¬ 
hibit an interesting property, leading to many curious 
effects: strongly spin-anisotropic interactions. Indeed, 
the interaction in the triplet configuration (same spins) 
is based on the exchange mechanism, and its strength is 
almost the same as for bare excitons, whereas in the sin¬ 
glet configuration an exchange leads to the formation of 
a dark exciton, whose energy is much higher than that 
of a lower polariton. It is therefore a second-order mech¬ 
anism, which is strongly suppressed. This concerns not 





9 


only interaction between two polaritons, but also between 
a polariton and a reservoir exciton. In both cases, the 
conclusion is that different potentials can be created for 


the two spin components using circular polarized optical 
pumping, either resonant (creating polaritons) or non¬ 
resonant (creating reservoir excitons). 


^ G. Volovik, The Universe in a Helium Droplet (Clarendon 
Press, Oxford, 2003). 

^ G. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

^ M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat 
Phys 2, 620 (2006), ISSN 1745-2473. 

^ O. Klein, Z. Phys. 53, 157 (1929). 

^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

® K. Novoselov, V. Fal’ko, L. Colombo, P. Gellert, 
M. Schwab, and K. Kim, Nature 490, 192 (2012). 

^ M. Katsnelson, G. Volovik, and M. Zubkov, Annals of 
Physics 336, 36 (2013). 

® G. Volovik and M. Zubkov, Annals of Physics 340, 352 
(2014). 

^ S. Zhou, G.-H. Gweon, A. Fedorov, P. First, W. de Heer, 

D. -H. Lee, F. Guinea, A. Castro Neto, and A. Lanzara, 
Nature Materials 6, 770 (2007). 

M. Polini, F. Guinea, M. Lewenstein, H. C. Manoharan, 
and V. Pellegrini, Nat Nano 8, 625 (2013), ISSN 1748-3387. 

L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu, and 

T. Esslinger, Nature 483, 302 (2012), ISSN 0028-0836. 

G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, 

T. Uehlinger, D. Greif, and T. Esslinger, Nature 515, 237 
(2014). 

P. Soltan-Panahi, J. Struck, P. Hauke, A. Bick, 
W. Plenkers, G. Meineke, C. Becker, P. Windpassinger, 

M. Lewenstein, and K. Sengstock, Nature Physics 7, 434 

( 2011 ). 

L. Duca, T. Li, M. Reitter, 1. Bloch, M. Schleier-Smith, 
and U. Schneider, Science 347, 288 (2015). 

O. Peleg, G. Bartal, B. Ereedman, O. Manela, M. Segev, 
and D. N. Christodoulides, Phys. Rev. Lett. 98, 103901 
(2007). 

U. Kuhl, S. Barkhofen, T. Tudorovskiy, H.-J. Stockmann, 
T. Hossain, L. de Eorges de Parny, and E. Mortessagne, 
Phys. Rev. B 82, 094308 (2010). 

E. Kalesaki, C. Delerue, C. Morais Smith, W. Beugeling, 
G. Allan, and D. Vanmaekelbergh, Phys. Rev. X 4, 011010 
(2014). 

Y. Plotnik, M. C. Rechtsman, D. Song, M. Heinrich, J. M. 
Zeuner, S. Nolle, Y. Lumer, N. Maikova, J. Xu, A. Szameit, 
et al., Nat Mater 13, 57 (2014), ISSN 1476-1122. 

A. B. Khanikaev, S. Hossein Mousavi, W.-K. Tse, M. Kar- 
garian, A. H. MacDonald, and G. Shvets, Nat Mater 12, 
233 (2013), ISSN 1476-1122. 

N. Goldman, 1. Satija, P. Nikolic, A. Bermudez, M. A. 
Martin-Delgado, M. Lewenstein, and 1. B. Spielman, Phys. 
Rev. Lett. 105, 255302 (2010), URL http://link.aps. 
org/doi/10.1103/PhysRevLett.105.255302 

M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, 
D. Podolsky, E. Dreisow, S. Nolle, M. Segev, and A. Sza¬ 
meit, Nature 496, 196 (2013), ISSN 0028-0836. 

D. Tanese, H. Elayac, D. Solnyshkov, A. Amo, A. Lemaitre, 

E. Galopin, R. Braive, P. Senellart, 1. Sagnes, G. Malpuech, 
et al.. Nature Comm. 4, 1749 (2013). 


D. Tanese, E. Gurevich, E. Baboux, T. Jacqmin, 
A. Lemaitre, E. Galopin, 1. Sagnes, A. Amo, J. Bloch, 
and E. Akkermans, Phys. Rev. Lett. 112, 146404 (2014), 
URL http://link.aps.org/doi/10.1103/PhysRevLett. 
112.146404 

T. Jacqmin, 1. Carusotto, 1. Sagnes, M. Abbarchi, D. Sol¬ 
nyshkov, D., G. Malpuech, E. Galopin, A. Lemaitre, 
J. Bloch, and A. Amo, Phys. Rev. Lett. 112, 116402 
(2014). 

A. Kavokin, J. Baumberg, G. Malpuech, and E. Laussy, 
Mieroeavities (Oxford University Press, 2011). 

C. Weisbuch, M. Nishioka, A. Ishikawa, and Y. Arakawa, 
Phys. Rev. Lett. 69, 3314 (1992), URL http://link.aps. 
org/doi/10.1103/PhysRevLett.69.3314 

J. Kasprzak, M. Richard, S. Kundermann, A. Baas, 
P. Jeambrun, J. M. J. Keeling, E. M. Marchetti, M. H. 
Szymanska, R. Andre, J. L. Staehli, et al.. Nature 443, 
409 (2006), ISSN 0028-0836. 

R. T. Whitlock and P. R. Zilsel, Phys. Rev. 131, 
2409 (1963), URL http://link.aps.org/doi/10.1103/ 
PhysRev.131.2409 

1. A. Shelykh, A. V. Kavokin, Y. G. Rubo, T. C. H. Liew, 
and G. Malpuech, Semiconductor Science and Technology 
25, 013001 (2010). 

N. Takemura, S. Trebaol, M. Wouters, M. Portella-Oberli, 
and B. Deveaud, Nature Physics 10, 500 (2014). 

A. Nalitov, G. Malpuech, H. Tercas, and D. Solnyshkov, 
Phys. Rev. Letters 114, 026803 (2015). 

V. G. Sala, D. D. Solnyshkov, 1. Carusotto, T. Jacqmin, 
A. Lemaitre, H. Tercas, A. Nalitov, M. Abbarchi, E. Ga¬ 
lopin, 1. Sagnes, et al., arXiv: 1406.4816 (2014). 

K. Bliokh, D. Smirnova, and E. Nori, Science 348, 1448 
(2015). 

C.-E. Bardyn, T. Karzig, G. Refael, and T. C. H. Liew, 
Phys. Rev. B 91, 161413 (2015), URL http://link.aps. 
org/doi/10.1103/PhysRevB.91.161413 
A. Nalitov, D. Solnyshkov, and G. Malpuech, Phys. Rev. 
Lett. 114, 116401 (2015). 

E. Wertz, L. Eerrier, D. Solnyshkov, R. Johne, D. San- 
vitto, A. Lemaitre, 1. Sagnes, R. Grousson, A. Kavokin, 
P. Senellart, et al.. Nature Physics 6, 860 (2010). 

C. Sturm, D. Tanese, H. Nguyen, H. Elayac, G. E., 
A. Lemaitre, 1. Sagnes, D. D. Solnyshkov, D. Amo, 

G. Malpuech, et al.. Nature Comm. 5, 3278 (2014). 

H. Nguyen, D. Vishnevsky, C. Sturm, D. Tanese, D. Sol¬ 
nyshkov, E. Galopin, A. Lemaitre, 1. Sagnes, A. Amo, 
G. Malpuech, et al., Phys. Rev. Lett. 110, 236601 (2013). 
C. Anton, T. C. H. Liew, G. Tosi, M. D. Martin, T. Gao, 
Z. Hatzopoulos, P. S. Eldridge, P. G. Savvidis, and L. Vina, 
Phys. Rev. B 88, 035313 (2013), URL http://link.aps. 
org/doi/10.1103/PhysRevB.88.035313 

A. Amo, T. Liew, C. Adrados, R. Houdre, E. Giacobino, 
A. Kavokin, and A. Bramati, Nature Photonics 4, 361 
( 2010 ). 

M.-H. Liu, J. Bundesmann, and K. Richter, Phys. Rev. B 
85, 085406 (2012), URL http://link.aps.org/doi/10. 



10 


1103/PhysRevB.85.085406 

D. -W. Zhang, Z.-Y. Xue, H. Yan, Z. D. Wang, and S.-L. 
Zhu, Phys. Rev. A 85, 013628 (2012), URL http://link, 
aps.org/doi/10.1103/PhysRevA.85.013628 

P. Wallace, Phys. Rev. 71, 622 (1947). 

N. Y. Kim, K. Kusudo, A. Lffler, S. Hfling, A. Forchel, 
and Y. Yamamoto, New Journal of Physics 15, 035032 
(2013), URL http://stacks . iop. org/1367-2630/15/i= 
3/a=035032 

E. A. Cerda-Mendez, D. N. Krizhanovskii, M. Wouters, 

R. Bradley, K. Biermann, K. Guda, R. Hey, P. V. San¬ 
tos, D. Sarkar, and M. S. Skolnick, Phys. Rev. Lett. 105, 
116402 (2010), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.105.116402 

N. Y. Kim, K. Kusudo, A. Lofiler, S. Hofling, A. Forchel, 
and Y. Yamamoto, Phys. Rev. B 89, 085306 (2014), 
URL http://link.aps.org/doi/10.1103/PhysRevB.89. 
085306 

B. Znahg, S. Brodbeck, Z. Wang, M. Kamp, C. Schneider, 

S. Hofling, and H. Deng, Appl. Phys. Lett. 106, 051104 
(2015). 

G. Panzarini, L. G. Andreani, A. Armitage, D. Baxter, 
M. S. Skolnick, V. N. Astratov, J. S. Roberts, A. V. Ka- 
vokin, M. R. Vladimirova, and M. A. Kaliteevski, Phys. 
Rev. B 59, 5082 (1999). 

A. Kavokin, G. Malpuech, and M. Glazov, Phys. Rev. Lett. 
95, 136601 (2005). 

G. Leyder, M. Romanelli, J. P. Karr, E. Giacobino, 

T. G. H. Liew, M. M. Glazov, A. V. Kavokin, G. Malpuech, 
and A. Bramati, Nat Phys 3, 628 (2007), ISSN 1745-2473. 
R. Hivet, H. Flayac, D. D. Solnyshkov, D. Tanese, 
T. Boulier, D. Andreoli, E. Giacobino, J. Bloch, A. Bra¬ 
mati, G. Malpuech, et ah, Nat Phys 8, 724 (2012), ISSN 


1745-2473. 

G. Dasbach, G. Diederichs, J. Tignon, G. Giuti, P. Rous- 
signol, G. Delalande, M. Bayer, and A. Forchel, Phys. Rev. 
B 71, 161308 (2005). 

V. Sala, D. Solnyshkov, 1. Garusotto, T. Jacqmin, 

A. Lemaitre, H. Tercas, A. Nalitov, M. Abbarchi, E. Ga- 
lopin, 1. Sagnes, et ah, Phys. Rev. X 5, 011034 (2015), URL 
http://link.aps.org/doi/10.1103/PhysRevX.5.011034 
P. Dirac, The Principles of Quantum Mechanics (Oxford, 
Glarendon Press, 1930). 

N. Hatano, R. Shirasaki, and H. Nakamura, Phys. Rev. A 
75, 032107 (2007). 

J.-S. Yang, X.-G. He, S.-H. Ghen, and G.-R. Ghang, Phys. 
Rev. B 78, 085312 (2008). 

H. Tercas, H. Flayac, D. Solnyshkov, D., and G. Malpuech, 
Phys. Rev. Lett. 112, 066402 (2014). 

B. Huard, J. A. Sulpizio, N. Slander, K. Todd, 

B. Yang, and D. Goldhaber-Gordon, Phys. Rev. Lett. 98, 
236803 (2007), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.98.236803 

A. Young and P. Kim, Nature Physics 5, 222 (2009). 

T. Salger, G. Grosser!, S. Kling, and M. Weitz, Phys. 
Rev. Lett. 107, 240401 (2011), URL http://link.aps. 
org/doi/10.1103/PhysRevLett.107.240401 
A. Galogeracos and N. Dombey, Int. J. Mod. Phys. A14, 
631 (1999). 

A. Galogeracos and N. Dombey, Gontemp. Phys. 40, 313 
(1999). 

C. Beenakker, Rev. Mod. Phys. 80, 1337 (2008). 

J. J. Pereira, F. Peeters, A. Ghaves, and G. Farias, Semi- 
cond. Sci. Technol. 25, 033002 (2010). 

P. Allain and J. Fuchs, Eur. Phys. J. B 83, 301 (2011). 


